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Abstract 

In this paper, we calculate the branching ratios for — )• D^tt^ — ?• Dfir^, — ?• D*^-k^ 
and B^ — )• D*~^7r^ decays in the perturbative QCD factorization approach. We find that the 
calculated branching ratios of these four decay channels agree well with the measured values and 
current experimental upper limit. In the numerical calculation, we take the decay constant and 
the shape parameter of the vector meson D* as fo* = 312 MeV and a/j* = 0.78 respectively, 
which are larger than those in the previous calculations. 
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I. INTRODUCTION 



In recent years, more and more effort has been made to the B meson decays with one 
[l| even two j2[ charmed mesons in the final states and it is found that the perturbative 
QCD factorization (pQCD) approach does work well in these decays. We will calculate 

the branching ratios for the B ~¥ Ds w decays, which are shown in figure 1, by employing 
the pQCD approach. The momenta of the two outgoing mesons are both approximately 
^7715(1 — {*)/^\)- This is still large enough to make a hard intermediate gluon in the 
hard part calculation. Most of the momenta come from the heavy b quark in quark level. 
The light quark u (d) inside meson, which is usually called spectator quark, 

carries small momentum of order of Aqcd- In order to form a fast moving light meson, 
the spectator quark need to connect the four-quark operator (bu)v-A{cs)v-A through 
an energetic gluon. The hard four-quark dynamic together with the spectator quark 
becomes six-quark effective interaction. Since six-quark interaction is hard dynamics, it 
is perturbatively calculable in theory. 

On the experimental side, the branching ratios of B^ — )■ Df7r~,B~^ — )■ -D^tt" and 
B^ — )■ D*~^7i~ have been measured by BaBar 13,] and Belle jl]. For — )■ D*~^7r^ decay 
only the experimental limit is given by CLEO 51]. We list their values in the following 

Br{B° D+n-) = (1.53 ± 0.35) x 10'^, 
Br{B^ ^ Djn^) = (1.6 ± 0.6) x 10~^ 
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Br{B^ ^ D*+7T-) = (3.0 ± 0.7) X 10^ 
Br{B^ Dl+n^) < 2.7 x 10^^ (1) 

This paper is organized as follows. In SectUll the light-cone wave functions of the initial 
and the final state mesons are discussed. In Sec lIIIt we calculate analytically the related 
Feynman diagrams and present the various decay amplitudes for the studied decay modes. 
The numerical results and the discussions are given in the section IIVI The conclusions 
are presented in the final part. 



II. WAVE FUNCTIONS OF INITIAL AND FINAL STATE MESONS 



In pQCD calculation, the light-cone wave functions are nonperturbative and not cal- 
culable, but they are universal and channel independent for all the hadronic decays. 

As a heavy meson, the B meson wave function is not well defined. In general, the B 
meson light-cone matrix element can be decomposed as 0] 



e'''^-^{0\UO)d^iz)\B{pB)) 



V2K 



(2) 



where n = (1,0, Ot), and v = (0, 1, Ot) are the unit vectors pointing to the plus and 
minus directions, respectively. Because the contribution of the second Lorentz structure 
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(psix.b) is numerically small and can be neglected, we only consider the contribution of 
the Lorentz structure: 

^b{x, b) = ^^(^B + meHM^, b). (3) 

In the heavy quark limit, we take the wave functions for the pseudoscalar meson Dg 
and the vector meson as 

^Dsix,b) = ^^75(^D. +mzjj0D,(x,6), (4) 
^Dtix, b) = i^D* + mDt)(pDt{x, b), (5) 

where the polar vector / = — (1, — '^'Ij^jOt)- In the considered decays, the D* meson 

is longitudinally polarized, so we only need to consider its wave function in longitudinal 
polarization. 

The wave function for the light pseudoscalar meson tt is given as 

$,(P, X, C) = ^^75 [J/'<Pni^) + ml<P^ {x) + Cm^(# - v ■ n)<Pl{x)] , (6) 

where P and x are the momentum and the momentum fraction of tt meson, respectively. 

The parameter ( is either +1 or —1 depending on the assignment of the momentum 
fraction x. The chiral scale parameter ttIq is defined as ttIq = m^/imu + m^). 



III. THE PERTURBATIVE QCD CALCULATION 

Using factorization theorem, we can separate the decay amplitude into soft, hard, and 
harder dynamics characterized by different scales, conceptually expressed as the convolu- 
tion. 



t) , (7) 



where A^j's are momenta of light anti-quarks included in each meson, and Tr denotes the 
trace over Dirac and color indices. C{t) is the Wilson coefficient which results from the 
radiative corrections at a short distance. In the above convolution, C{t) includes the 
harder dynamics at a larger scale than that at the Mb scale and describes the evolution 
of local 4- Fermi operators from mw (the W boson mass) down to t ~ scale, 
where A = Mb — ttli,. The function H{ki,k2,k3,t) describes the four quark operator 
and the spectator quark connected by a hard gluon whose is in the order of AMb, and 
includes the 0{\/AMb) hard dynamics. Therefore, this hard part H can be perturbatively 
calculated. The function ^ the wave functions of Dg and tt. 

Since the b quark is rather heavy, we consider the B meson at rest for simplicity. It 
is convenient to use the fight-cone coordinate {p^,p~ , Pr) is used to describe the meson's 
momenta: 

p± = i=(/±p3)^ and PT = (p\/). (8) 
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FIG. 1: Diagrams contributing to the decays B ^ Dg it. 

At the rest frame of B meson, the hght meson moves very fast and so or can 
be treated as zero. Using these coordinates, the B meson and the two final state meson 
momenta can be written as 

Mr Mr o Mr 

P5 = -|(1,1,0t), P2 = ^(l,r^0T), P3 = ^(0,l-r2,OT), (9) 

respectively, where r = M (*) /Mb- Putting the hght anti-quark momenta in Dl and 
71 mesons as ki, k2, and fcs, respectively, we can choose 

fci = (xiP+,0,kiT), k2 = {x2P^,0,k2T), A;3 = (0,X3P3-,k3T). (10) 

For these considered decay channels, the integration over k^ , /c^, and k^ in equation ([7]) 
will lead to 



dxidx2dx'ihidhih2dh2hzdhz 



•Tr 



where 6j is the conjugate space coordinate of kiT, and t is the largest energy scale in the 
function H{xi,bi,t). The last term e"'^^*-' in equation f|TT]) is the Sudakov form factor 
which suppresses the soft dynamics effectively [8|. 

For the considered decays, the related weak effective Hamiltonian H^ff can be written 
as 

n^ff = ^ mif^Wf,) + C2{f^)02{fi))] , (12) 

where the four-quark operators are 

Oi = {baUfj)v-A{CaS^)v-A, O2 = {baUa)v-A{CaSa)v-A, (13) 

with a, (3 being the color indexes, and {qiq2)v-A = Qil^i^ — 1^)<12- The Fermi constant 
Gp = 1.16639 X 10~^GeV^~^ and Ci,2(At) are Wilson coefficients running with the renor- 
malization scale /i . The leading order diagrams contributing to the decays B — > Pi*V 
are drawn in figure [T] according to this effective Hamiltonian. 
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In the following, we will get the analytic formulas by calculating the hard part H{t) 
at leading order. Involving the meson wave functions, the amplitude for the factorizable 
tree emission diagrams Fig. 1(a) and (b) can be written as: 



pL poo 

Fe = SirCpf M / dxidx3 / feirffoi 636/63 ^^(xi, 61) 
^ Jo Jo 

X { [(0:3 + 1)0^(X3) - r.(2x3 - iMix,) + 0^x3))] 
xEe{t)he{xi, 3:3(1 - bi, h)St{x3) 

+2rX{xs)E,{t')K{x,,x,{l - rl,,,),b,,h)St{x,)} , (14) 

where Cp = 4/3 is the group factor of SU{3)c gauge group, and the mass ratios = 
fn^/mB^r (*) = m (*)/mB- Here / (,) is the decay constant of d'^^ meson, and St{x) is 
the jet function [10]. The factor evolving with the scale t is given by: 

E,{t) = as{t) exp[-SB{t) - S^{t)], (15) 



where S'B(t), S'^(t) are expressions for Sudakov form factors [10||. The hard function is 
written as 

he{Xl,X2, bi, 62) = Kq {y/x^niBbi) [9{bi - b2)Ko (y/x^niBbi) Iq {y/x^mBb2) 

+9{b2 - bi)Ko {y/x^mBb2) Iq {^/x^mBbl)] . (16) 

The hard scales t^'-* in Eq.f ll4p are determined by 

t = max{^x-i{l - r'^^i^,^)mB, I/61, I/63), 

t' = max{^xi{l - r'^^(,))mB, I/61, 1/63)- (17) 

For the nonfactorizable tree emission diagrams Fig. 1(c) and (d), all three meson wave 
functions are involved. The integraton of 63 can be performed using 5 function 5(63 — 62) 
and the result is 

pl poo 

Me = ~16n^/2NcCF / dxidx2dxs / 6i(i6i 62C^^2 $5(2^1, ^i)^^{*) (2^2) 

Jo Jo 

X { [{x2 - l)(p^{x3) + r^X3(0f (X3) - (plixs))] En{t)h\{xi, X2, 0:3, 61, 62) 

+ [(x3 + X2)0;^(x3) -r^X3(0^(x3) + 0^(x3))E„(t')/i^(a;i, X2, X3, 61, 63)] } , (18) 

where the expressions for the evolution factor is En = as(t) exp[—S(t)\b^=b-^^] with the 
Sudakov exponent S = Sb + + S.,,. 

The hard functions /i^, i = 1, 2 in the amplitude are given as 

hi, = [9ib^-b2)Ko{Ab^)Io{Ab2) + 9ib2-b^)Ko{Ab2)Io{Ab,)] 



Ko(Gi62), for G'2>0 



^ifMVmbs), for G^<0, . . 
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TABLE I: Input parameters used in the numerical calculation [q. 1 1 1 1 . 



Masses 


= 0.14 GeV, = 1.3 GeV, 
niD, = 1.9685 GeV, mo* = 2.1123 GeV, 
mB = 5.28 GeV, raw = 80.4 GeV, 


Decay constants 


/b = 0.19 GeV, U = 0.13 GeV, 
/d. = 0.273 GeV, fnt = 0.312 GeV, 


Lifetimes 


Tb± = 1.638 X 10-12 ^ i_53Q iq-12 


CRM 


Vcb = 0.0412 ± 0.0011, Vus = 0.2255 ± 0.0019. 



with the variables 



A' 



XiX3(l - r^(,))m|, 



= (xi + X2)r^(.) - (1 - Xi - X2)x3(l - r'^^(,))m 



By 



Gl 



{xi -X2)x3(l -r^(,))m|. 



The hard scales in Eq.f lT^ are given by 



t = max(AmB, y^Gfrns, 1/62), 
t' = max(ylmB, -y/G^ms, 1/61, 1/62). 



Then the total decay amplitude oi B ^ Ds n decays can be written as 



(20) 



(21) 



(22) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



For the numerical calculation, we list the input parameters in Table HI 
For the B meson wave function, we adopt the model 



(f)B{x,b) = Nbx'^{1 — x^ex-p 



Ml x2 1 



(23) 



where ujh is a free parameter and we take ujb = 0.4 ± 0.04 GeV in numerical calculations, 
and Nb = 91.745 is the normalization factor for Ub = 0.4. 
For oi*^ meson, the distribution amplitude is taken as: 



!,^(,)(a;) = _ x) 



1 



1 - a^(,)(l - 2x) 



(24) 



with the Gegenbauer coefficients ao^ = 0.3 and qd* = 0.78. The CLEO and BarBar 
collaborations reported their work on the measurements of the decay constant of Dg 
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meson and obtained fo, = 274 ± 13 ± 7 MeV [lH and 283 ± 17 ± 7 ± 14 MeV [jj, 
respectively. However, the decay constant of the vector meson has not been directly 
measured in experiments so far. From the conclusions draw by the CLEO collaboration 
[111 , one can find that there exists a relation: 



fn* Id Ib 



[1.1,1.2], 



(25) 



which is consistent with that from lattice simulation [13| and the QCD sum rules calcula- 
tions IJ]. From ta ble p it is easy to see the value of the ratio fn* / fn^ is 1-14 in our work. 
It is different from [l5|, where the relation between fn* and fn, derived from HQET 



fOs 



rrin: 



(26) 



was used. From this equation, one can get the value of fn* , which is smaller than that of 
fos- 

The twist-2 pion distribution amplitude (p^, and the twist-3 ones (p^ and 0^ have been 
parametrized as 



f. 



zQxil — x) 



2V2iVc 
+alCl'\2x-l) 

f. 



1 + alCl^\2x - 1) + 4^2^" (2a; - 1 



3/2, 



1 + (30r73 - ^pi)Cl^\2x - 1) - 3 <; rj^us + --p^(l + 6^ 



2v/2iVe 
xCy^(2x- 1) 

:(l-2a;) 
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20' 



2V2iV, 



1 

1 + 6(5?73 - -rjsus 



pX)(1 - 10a; + lOx^) 



(27) 

(28) 
(29) 



with the mass ratio p,r 
Cf/2(t) 



2C 5' 

rriu + md)/m.n = m.^/mQ and the Gegenbauer polynomials 



1 
2 

3t 



(3^2 - 1) 



Cy\t) = -(3-30t' + 35t^), 



1 



(t) = ^(5t2-l). 



— (l-14t2 + 21t^). 



The Gegenbauer coefficients are given as 



= 0, 4 = 0.115, al 



-0.015. 



(30) 
(31) 
(32) 

(33) 



The values of other parameters are taken as [l6j 773 = 0.015 and u = —3.0. 
In the B-rest frame, the decay width of B ^ Dg vr can be obtained by 



327r 



(34) 



TABLE II: Branching ratios (xlO~^) for the decays Dfir^ , D*^tt^ and i?+ 

Df'K^,D*^'ir^. The first theoretical error is from the the B meson shape parameter a;^. The 
second error is from the higher order pQCD correction. The third one is from the uncertainties 
of CKM matrix elements. 



Channel 


This work 


Data 


D+-K- 
B+ ^ D+vrO 
B^ D*+-K- 
B+ D*+7r^ 


1 ot;+0. 36+0. 41+0. 10 
■•-•"^-0. 52-0. 56-0. 10 
1 QQ+0. 39+0. 81+0. 11 
-•-•^"-0. 56-0. 31-0. 11 
9 cQ+0. 45+0. 70+0. 15 
^•"-'^-0.76-0.60-0.15 
9 7O+0. 48+0. 74+0. 16 
-^^ 'O-0. 82-0. 65-0. 16 


1.53 ±0.35 
1.6 ±0.6 
3.0 ±0.7 

< 27 




FIG. 2: Branching ratios (in units of 10 ^) of 5° D*^tt and B^ D*+7r° decays as 
functions of Gegenbauer moment a^* ■ 



where A is the total decay amplitude shown in Eq. fl22|) . 

Using the wave functions as specified in the previous section and the input parameters 
listed in this section, it is straightforward to calculate the CP-averaged branching ratios 
for the considered decays, which are listed in Table [TTl The first error in these entries is 
caused by the B meson shape parameter Ub = 0.40 ± 0.04. The second error is from the 
higher order pQCD correction: the choice of hard scales, defined in Eq.f llTp and Eq.( pT]) . 
which vary from 0.9t to l.lt. The third error is from the uncertainties of the CKM matrix 
elements which are listed in table [H 

In previous calculations jH ^ , the authors have considered that the value of the Gegen- 
bauer moment a£,* was the same as that of a^^ and taken them as 0.3. Here we take 
O'D* = 0.78, which is determined to fit the requirement that (f)D*{x), shown in Eq. f l24p . 
has a maximum at x = . In Fig. [21 we plot that a^,* dependence of the branching 

ratios of — )■ D+*7r~ and — )■ D^*7r'^. One can find that the branching ratios are not 
sensitive to the variations of qd*- 

From the numerical results, we find that the non-factorizable contributions are very 
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small and almost neglectable. They are about 10% of the factorizable ones in each decays. 
The main contributions come from the factorizable amplitudes. 



V. CONCLUSION 



In this paper, we calculate the branching ratios of decays — t- D^ti , — t- D'^tt^, 
B^ — )> -D*+7r" and 5+ — > D*+7r° in the pQCD factorization approach. We find that: 

• The decays considered here have branching ratios about 10^ smaller than those of 
the B — > D^*^7r decays, and they comes mainly from the relevant CKM matrix 
elements. 

• From the numerical results shown in table[Tll one can find that the pQCD predictions 
for these considered decay channels are consistent with the measured values and 
currently available experimental upper limit. 

• To determine decay constant of the vector meson D*"*" , the relation 

Id* fo fs 

is used. It indicates that the value of fo* is larger than that of fn* , which is contrary 
to the conclusion derived from the relation 



In the numerical calculation, we take an* = 0.78, which is larger than the value 
given in the previous calculations. It is determined to fit the requirement that the 
wave function (f)£)*{x) has a maximum at x — "'Os-mc 
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